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ABSTRACT 


An improved model for a once-through counter-cross-flow 
waste heat recovery unit with a segmented fin-tube arrange- 
ment was developed. The model was coded in FORTRAN IV 
computer language for use on an IBM-3033 computer system, 
and was tested for various conditions of uniform and non- 
uniform gas glow distributions. Additional parameters 
varied were heat exchanger size, steam outlet pressure, 
steam outlet temperature and gas flow rate. The effect of 
each variation on WHRU performance was evaluated and the 
interrelationship discussed. Results indicate the ideal 
design is a high steam temperature, low operating pressure, 


flow distribution controlled heat exchanger. 
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tee LN eRe DUCTION 


A. BACKGROUND 

mmcceeme IMtroductionwef the DD-968 class ships into 
the fleet, the U.S. Navy has become firmly committed to gas 
turbine propulsion plants. Currently, besides the DD-963 
class, there are four additional ship classes which are, or 
will be, gas turbine powered. These are the Patrol Hydro- 
foil Missile (PHM's), the Oliver Hazard Perry FFG-7 class, 
the Ticonderoga CG-47 class and the soon to be built DDG-51 
class. 

This commitment to gas turbine propulsion has stimulated 
great interest in the concept of converting the waste heat 
from the gas turbine engines into useful shaft power. In 
December 1979 the Naval Sea Systems Command awarded contracts 
for development of a conceptual design of a new generation 
of waste heat recovery systems. The goal of the design was 
to provide a system which would realize a 25 percent minimum 
improvement in ship propulsion fuel consumption and increase 
the ship's operating range by 1000 nautical miles. The out- 
growth of the effort has been the RAnkine Cycle Energy Recovery 
(RACER) system. As conceived, the RACER system will be an 
unfired waste heat recovery system designed to convert waste 
heat from the exhaust of the main propulsion gas turbine 


engines into useful shaft horsepower. Such a system has 


Ais, 





commonly been called a COGAS or COmbined Gas And Steam system. 
The hot exhaust gases of the main propulsion turbine are 
passed through a Waste Heat Recovery Unit (WHRU), where steam 
of the desired temperature and pressure 1S generated. The 
steam is used to drive a steam turbine which operates in 
parallel with the propulsion gas turbine through a common 
reduction gear. 

Numerous papers have been written concerning the RACER 
system. Halkola, Cambell and Jung [{Ref. 1] discussed the 
conceptual design for the system. Mattson [Ref. 2] presented 
a plan for designing reliability and maintainability into 
the RACER system and Combs [{Ref. 3] developed a computer model 
of a COGAS system which designs a waste heat recovery unit 
and allows comparison of system performance by variation of 
system parameters. 

References 1, 2 and 3 all agree that the desired goal of 
fuel conservation must not come as a result of increased 
Maintenance, manpower, or main propulsion system down-time. 
The design philosophy dictates that the system have high 
reliability, be inherently safe, have low maintenance require- 
ments and be essentially automatic in operation. 

Reference 1 states that the waste heat recovery unit 
paces the design of the system. The majority of effort by 
Combs [Ref. 3] was to create a WHRU design model which would 
Simulate the performance characteristics of a physical boiler. 


Since great expense is involved in the design, fabrication 
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and operation of prototype systems, it is of primary impor- 
tance to achieve a realistic mathematical simulation first. 
To this end, Combs devised a model which designs a WHRU 
when given feedwater inlet temperature, steam outlet tempera- 
ture, operating pressure, inlet and outlet gas conditions 
and gas flow rate. 

In this study, the initial model of Combs was considered 
and an improved computer model of the waste heat recovery 
unit was developed. A once-through counter-cross-flow heat 
exchanger was selected because of the advantages of lower 
initial cost, faster response, compactness and low operating 
cost. The improvement in the model is that the heat exchanger 
is discretized into individual "nodes" much in the manner 
done by Shu and Kuo [Ref. 4]. This procedure treats smaller 
physical pieces of the overall heat exchanger as individual 
small heat exchangers to attain a better understanding of 
the temperature distributions within the WHRU. The improved 
model allows an analysis of the effect of non-uniform velocity 
distributions on the performance of the WHRU, and ultimately, 


the RACER system. 


Bee OBJECTIVES 
There were five main objectives for this thesis. 
1. Devise an improved computer model for the once-through 
WHRU. This model was to predict the performance of various 


size WHRU's given feedwater inlet temperature, steam outlet 
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temperature, WHRU operating pressure, inlet gas mass flow 
distribution and inlet gas temperature distribution. 

2. Enable the designer to establish a design WHRU size 
based on the performance of the model with various input 
Senditions. 

3. Evaluate the design WHRU utilizing various uniform 
and non-uniform inlet gas distributions for both on-design 
and off-design conditions. 

4, Re-evaluate the WHRU utilizing the same uniform and 
non-uniform gas distributions into the WHRU, but with 
different steam outlet temperatures. 

5. Provide a framework for future studies involving a 
development of a three dimensional model of a WHRU and a 


detailed, complete COGAS system cycle analysis and optimization. 
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A. OVERVIEW 

The waste heat recovery unit chosen was a once-through 
counter-cross-flow heat exchanger. A computer program was 
written to model this type heat exchanger for analysis of 
the effect of uniform and non-uniform inlet gas mass distri- 
butions. Variation of other parameters, other than the gas 
distribution, provides a thorough understanding of the oper- 
ation of the heat exchanger in a multitude of conditions. 
The model was written for incorporation into a RACER system 
analysis, but it is not limited to that application. 

Figure 1 shows a three dimensional representation of a 
typical layout of a marine gas turbine waste heat recovery 
unit. The unit is headered on the fluid inlet and exit and 
designed to be installed in the gas turbine exhaust gas 
duct. The exhaust gas entering the bottom of the WHRU pro- 
vides the energy for transfer to the fluid-side of the system. 
The fluid (feedwater) enters at the upper header of the WHRU 
and flows alternately back and forth through the unit until 
exiting as steam at the lower header. The U-bends shown in 
Figure 1 at the ends of the tube passes allow each tube to 
enter and exit the heat exchanger as a continuous path for 
fluid flow. As a tube makes one horizontal traverse across 


the WHRU it is referred to as a pass. Therefore, the number 
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of times a tube passes through the gas flow defines the 
vertical size of the heat exchanger as a 15, 20, 25, etc., 
pass heat exchanger. 

A gas flow distribution is depicted as entering the WHRU 
in Figure 1. The gas flow is shown oriented along the tube 
pass aS a non-uniform distribution. Additionally, the gas 
flow may be non-uniform in the direction into the heat 
exchanger, 1i.e., along the axis of the headers as shown in 
Figure 1. This two dimensional gas flow distribution may be 
non-uniform in both mass and temperature. The nature of the 
heat exchanger, by virtue of the cross flow of the fluid, 
decrees that, even if the mass and temperature distributions 
entering the heat exchanger were uniform in both directions, 
these distributions can become distorted and non-uniform as 
the gas passes through the heat exchanger in the third 
dimension. This distortion can be caused by the variation 
of the heat transfer characteristics along the tube. The 
problem of analysis becomes more complicated if cross mixing 
of the gas is considered. 

The model of Combs [Ref. 3] treated the gas distribution 
as a uniform and constant value in one direction. Combs 
assumed the mass and temperature of the gas along the axis 
of the tube was everywhere the same for a complete pass. 

The model yielded results which were close to the results 
attained by other studies and gave a reasonable estimate of 


the overall heat transfer characteristics of the heat 
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exchanger. The model, however, could not predict the effect 
of non-uniformity of the gas distribution on the performance 
of the WHRU. This study discretizes the dimension along the 
axis of the tube pass and enables the designer to evaluate 
the effects of variations in mass flow rate and temperature. 
It 1s assumed, however, that the distributions along the 
axis of the headers are all the same as the outer distribu- 
tion. The concept of discretization and "nodal" arrangement 
is shown in Figure 2. The representation shows a matrix 
arrangement where the number of rows corresponds to the 
number of passes in the heat exchanger, and the number of 
columns corresponds to the number of gas paths in the 
discretization. The number of passes is a user defined input 
to the model but the number of gas paths is fixed in the 
model at ten. 

Use of the model requires that the user provide the 
fPeeevowlng conditions. First, the number of passes that a 
particular water/steam tube is to make through the heat 
exchanger must be specified. This parameter will physically 
Size the heat exchanger in the vertical direction. The mass 
flow rate of the gas turbine exhaust gas, (m,) , for each of 
ten exhaust gas paths through the heat exchanger must be 
entered. For a uniform distribution, this is simply the 
total exhaust gas mass flow divided by ten. A simplifying 
assumption of the model is that there is no cross mixing of 


flow as the gas passes through the heat exchanger; therefore, 


fa 





the mass entering a particular gas path is constant through 
the heat exchanger. 

Next, the temperature of the gas entering each gas path 
is entered. For this study the temperature distribution of 
the gas waS considered uniform for all runs; however, the 
model is not constrained to uniform temperature distributions. 

The gas was assumed to be at atmospheric pressure 
entering the heat exchanger and gas side pressure drops were 
not considered in the calculations. 

The user next enters the feedwater inlet temperature, 
the steam outlet temperature and the fluid-side operating 
Peescscure. These parameters are used in the initial calicula- 
tions to determine a first "guess" mass flow rate on the 
fluid-side of the heat exchanger. The minimum allowed 
exhaust gas temperature is entered. This value 1s a user 
desired constraint and is specified as the minimum acceptable 
temperature the user can accept to prohibit hot side corro- 
Sion. The model is not constrained to maintain the exhaust 
gas temperature above the user specified temperature. 

The number of feet in length, per node, is entered. 

This parameter is based on the physical size of the heat 
Semeanger being analyzed in the horizontal direction. For 
Pecans Im this study, an assumed tube row length of ten 
feet divided by ten gas paths yielded a tube length per node 
of one foot. There is no constraint on the length parameter, 


within reasonable limits. 
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The heat exchanger model was written as a two dimensional 
analysis, 1.e., all parameters in the third dimension were 
assumed uniform. The next parameter entered is the number 
of tubes per row, which is used to size the heat exchanger 
in the third dimension. Figure 1 shows a three dimensional 
representation of the heat exchanger design. Each particular 
tube enters at the top of the heat exchanger and passes 
alternately back and forth below itself until exiting the 
heat exchanger at the bottom. The tubes in each subsequent 
plane, behind the first, enter and exit in the same manner, 
but alternate slightly up and down to give a staggered tube 
arrangement. The nodes are numbered in standard matrix format 
with the upper left being node (1,1), and counting across 
Mmevemone as (1,2), (1,3), etc., to (1,10); then node (2,1) 
is directly below node (1,1), and so on, to the bottom of 
the matrix. Figure 2 shows the nodal arrangement of the 
heat exchanger. The fluid may enter from either side of 
the heat exchanger. If the heat exchanger has an even number 
of passes, the fluid will enter on the same side as the 
steam exit, and on the opposite side for an odd number of 
passes. Each tube is continuous from entry to exit. 

The final user entry is the geometric scaling factor 
which is used in the geometry portion of the model to scale 
all the physical dimensions of the tube, with the exception 


@eeeene length, up or down. 
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Two additional simplifications were assumed in the model. 
First, the fluid side pressure drop was not considered. 
Lastly, even though the tube passes were discretized to 
determine the temperature distributions, the tube length 
within a node still has a significant length for a tempera- 
ture gradient. For analysis, the gas temperature at any 
node was assumed to be uniform across the length of the node 


at both the entry and exit. 


Pee LNITLAL CALCULATIONS 

The WHRU is essentially analyzed in three sections; 
Superheating, boiling and heating. Given an inlet gas distri- 
bution, water/steam inlet and outlet temperatures, and an 
operating pressure, there is a unique fluid-side mass flow 
rate which will satisfy all boundary conditions. The model 
eeeemets, through iteration, to determine the resultant 
fluid-side mass flow rate. The analysis is also greatly 
affected by the number of passes (i.e., the physical size) 
of the heat exchanger. There are clearly combinations of 
parameters for which no solution is possible, or the solu- 
tion is trivial in that the mass flow rate of the fluid is 
so small as to be meaningless. 


For calculations, variable assignment is as follows: 


(1) m_{m,n) = gas flow rate through gas path n at tube 
J pass ™; 

eel ed) = Gas temperature into node (m,n). The 
3 gas temperature out of node (m,n) will 


be the gas temperature into node (m-l,n); 
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(3) T, (m,n) = fluid temperature into node (m,n). The 
fluid temperature out node (m,n) will 
equal the fluid temperature into node 
(m,nt+l) if the fluid is flowing to the 
migdesee fig. 3), Or the fluid 
temperature out of node (m,n) will equal 
the fluid inlet temperature of node 
Gio ete ethe fiuid is flowing to the 
left. 

At the end of a fluid pass the fluid 
out of a node will equal the fluid 
temperature into the node directly below, 


on T, (m,n) = T.(mtl,n); 


(4) M = fluid-side mass flow rate 

Figures 3 and 4 are representations of typical nodes 
depicting the above explanation of fluid and gas flow. 

After all entries are made the program initializes by 
summing the individual m,(NR,n) elements to find mM sotar’ 
then the ey elements are averaged. The Be eae ele- 
ments are averaged by multiplying each element by the corres- 
Pemeing fraction of the total gas mass that is passing through 
each gas path. These "weighted" temperatures are summed and 
the result is divided by the number of gas paths. The result 
is then the average temperature of the gas into the heat 
exchanger. An average gas temperature is determined using 


the user specified minimum gas temperature out and the aver- 


age gas temperature into the heat exchanger. Then: 


_ - ae wotd 1 Cag G avgin — tT avg out?’ 


where Cog is evaluated at the heat exchanger gas bulk 
temperature. Now, a first "guess" fluid mass flow rate is 


Calculated from an energy balance 


Ie) 





Q 
h 


f OE oats oa yi 


where rs and re are the fluid inlet and outlet 
enthalpies at the user specified feedwater inlet temperature 


and steam outlet temperature. 


ieee GEOMETRY 

The geometry portion of the model is nearly identical to 
that used by Combs [Ref. 3]. A fin-tube with helically 
curved extended surfaces on Circular tubes is assumed for 
this model. The fins are segmented and the tubes are 
configured in banks of one row each with the rows staggered. 


The base tube surface is taken from Weirman, Taborek and 


Marner [Ref. 5]. The description of the finned tubes is as 
follows: 
d. = tube inside diameter = ae 86. - in. 
qd = tube outside diameter eee OO. 
ds = fin root diameter eee Lie 
Ne = fins per inch = 5.94 
Q = fin height =e sty, 
dn Sere matn Of CUE £rom fin tip = Oe 82 in. 
d = fin outside diameter = 4,03 in. 
te = fin thickness = 0.048 in. 
= fin segment width =~ eo. 17 sin. 
ae = number of segments in 360 degrees = 38 


26 








The finned tube configuration 1s shown in Figure 5. The 


tube length, L, and the number of tubes per row, N are 


1 a aid 
the user specified entries. The tube layout is shown in 
Figure 6. The center to center tube spacing in the trans- 
verse direction is 4.50 inches. The equilateral, staggered 
tube arrangement used in this model leads to a spacing 
normal to the gas flow, Su: of 4.50 inches and a spacing 
parallel to the gas flow, sei of 3.90 inches, as shown in 
Figure 6. The heat exchanger height is established by the 
number of WHRU passes (rows) and the tube spacing parallel 
to the gas flow, aoe 

In order to establish the minimum gas flow cross-sectional 


area the total "blocked" frontal area, Avs of the heat 


exchanger must be calculated from 


A = N 


iG + iG NN ae a P 
b Oo 


E/ xr je a go = 


and the minimum gas flow area is 


where A. is the total frontal area of the heat exchanger 
exposed to gas flow. The total inside area available for 


heat transfer per pass is 


Pe = i dq. iB Nest 


2), 





The gas-side surface area available per pass for heat 
transfer is the sum of the fin surface area and the bare 
tube surface area per pass. The fin surface area per tube 


is calculated from 


Z ie. = aa) 
Aci, = N,LINS(2 2, wl 5+ 2t, 2 + wit.) + s(da, - a ] 
where dey - di. = 2h. The bare tube area per tube is 

Aye = Td, hb- td t, Nye b 


Therefor2, the total outside area per pass available for 


heat transfer is 


Bop 7 Mie poyleyas i Ane? 
The cross sectional area for fluid flow is calculated from 


2 


qs Ne/p 


D> 
| 4 


ia 


With the mass flow rates, terminal temperatures and heat 
exchanger geometry established, the remainder of the model 
may be solved by nodal analysis for the unique fluid mass 
E¥ew rate and the location of the beginning of the boiling 


and superheating regimes. 
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D. SUPERHEATER NODAL ANALYSIS 

In the heat exchanger there is only one node where two 
conditions are initially known, and that is at the node 
where the steam is exiting. At this node, both the gas inlet 
mass flow rate and temperature and the fluid outlet tempera- 
ture have been specified by the user. Using the initial 
values of the fluid mass flow rate from the estimated overall 
energy balance and the desired superheater steam outlet 
temperature, the heat exchanger is analyzed in the reverse 
direction through the superheater, through the boiler and 
through the heater to arrive at the feedwater inlet. For 
the given geometry and flow conditions, the model will calcu- 
late the water inlet temperature consistent with the assumed 
estimated mass flow rate. If the temperature arrived at by 
analysis as the inlet to the first node does not match the 
user specified feedwater inlet temperature, then the fluid- 
Side mass flow rate is adjusted up or down and the entire 
process is begun again. 

To begin, at the final node the gas inlet temperature 
and mass flow rate are known. Also, the fluid outlet tempera- 
ture and a first "guess" fluid mass flow rate are known. 
It is then assumed that the bulk (average) fluid and gas 
temperatures for the node are equal to these known tempera- 
tures. All fluid and gas properties for a node are evaluated 
at the bulk temperatures. 

The gas-side Reynolds number for the node is calculated 


uSing the bulk temperature. With the Reynolds number, a 
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J=Lactor is obtained from a polynomial fit to the data for 
mepe layout number 5 in Ref. 5 (Fig. 6). By definition, 


the j-factor is related to the heat transfer coefficient, 


oo by 
a 2/3 
Pye te roducing 
_ Nu 
2S) = Semmes 


the previous expression can be written as 


oF Nu pr-/? _ Nu 
J Re Pr ~ Wye’ 
g g Re Pr 
g 
and 
Nu = j Re pri/3 
g g 
where 
le el 
Nu = 3 
g 


and Be SeeicmrEMermal conductivity of the gas. Therefore, a 
Fefationship may be written for the heat transfer coefficient 


as follows: 
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= 

iI 

ee 
Dee 


Re Pr 
g 


In the superheat region, the water-side heat transfer 
coefficient is calculated by using the Dittus-Boelter 


correlation 


where all properties are obtained at the bulk temperature of 
the steam in the node. Using the tube wall resistance 


together with Ae and h the overall heat transfer coeffi- 


£’ 


cient for the node can be written as 


i 








U5j (men) = | wt San SEY A. ‘ 
h i a7 K N : i a S h 
f We sD op Ne g 
where: 
A.. 
- - fea 
Op 


iiesiin efficiency, Nee is calculated from the expression 


tanh ML 
Ne = ML 
where 
i a vha P/KA 


Sell 





fod, e2t eae fin 1s approximated by a set of rectangular 


strips extending from the tube wall, 


Now, the cross-sectional area, A, of a rectangular strip 


fey be written as 


and the perimeter is 


eemmusing Known geometric parameters and the thermal conduc- 


tivity of the fin metal, the parameter ML may be expressed 


as 
[ie — Ou 1) 
where 
ao = -~ d , 2(wott,) 
2 Wo He ae 
S SE 
Now, 
tanh Cvh Acsn 
Dees t+ = (1 - — ) ] 





The pass effectiveness is calculated using the effective- 
ness-NTU method. The correlation for the effectiveness for 
a single pass, cross flow heat exchanger with both fluids 


unmixed is taken from Incropera and DeWitt [Ref. 6] 


ee = J- exp [ (a) (NTU) °*** {exp {-c_(wru) °° 7%) -13 
r 
where: 
Cmin 
oa ae 
max 


NTU = — ~~ 
‘enm 
min 
co. 1s the lesser of C. and C_ where 
min if g 
Ce : ist ne 
and 
C = C_ m_(m,n) 
g PG J 


and ors, 1s the greater of the two. All fluid and gas 
properties are again evaluated at the bulk temperature. 
Two equations must be satisfied to find the solution to 


a node. These are the rate equation 


hea C11 12 okie. 
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and the energy equation 


2 = Me Saye ae out Te TA 


In the two equations the only unknown is T vege Equating 
the two and solving for T. oe yields 

< _ Sani Ty ne we C of Te out 

1c gl SC t : 


aim | es Si 


Now, Q for the node is 


e “pete out ~ Te in 


and 


ab = T_. _ 
G-Oue gin 


Ty bulk — gin Ty out 


and 
“eee 2 ie oes gayle 


An average cutside tube wall temperature is calculated from 
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Cae ; 
T = T eee Ely 
WO g bulk De A 
where A, /ALS are total inside/outside areas for the node 
and toe? Mae are the new node bulk temperatures. The expres- 
Sion for ee is derived from the following formulation for 


heat transfer in the node 





R 
— = O = 
oS To yR,, gB T ep) 
where 
lL 
QR = 
Oo Ny Ne Ay, 
and 
_ 1 
Desig 
or “tz 


: 2 Se tie 
25 2 
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This gas-side film temperature is now introduced at the 
beginning of the calculations for the node as the new gas 
bulk temperature. 

The assumed bulk temperatures, at the beginning of the 
calculations, were obviously too high since the gas inlet 
and fluid outlet temperatures were used. This assumption, 
however, allowed a solution for the node to be found. The 
model stores the result of this first pass and returns to 
the beginning utilizing the new bulk temperatures to calcu- 
late a new solution to the node. At the end of the second 
is 


pass through, the result of the calculation for terete 


Gempared to the result obtained the first time. If Ee wire 
has changed, then the new result is stored and the process 

1s repeated. When We Gut: no longer changes, the process has 
converged to a good solution for the node. 

A check is made to determine if the fluid inlet tempera- 
ture is less than or equal to the saturation temperature for 
the pressure specified. If the temperature is higher than 
the saturation temperature, the model proceeds to the next 
node of the superheater and repeats the entire process. 

If the fluid temperature equals the saturation temperature, 
Superheating has theoretically begun at the interface between 
Enis and the final node of the boiling portion of the model. 


If this condition exists, the model proceeds to the boiling 


nodal analysis. 
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As mentioned previously, the model works through the 
feat exchanger in the direction opposite to fluid flow. 

When it is determined that the fluid temperature into a node 
is below the saturation temperature, the node is "Split" to 
determine the physical location in the node where super- 
heating begins. 

To split the node the m quantity into the node and the 
physical length, L, are halved. The geometry is recomputed 
as if the node were only half as large. The solution 
process is reinitiated and, when a solution is achieved, the 
Saturation temperature check is redone. If the temperature 
into the node is still below the saturation temperature, m 
and the length are halved again and the process is repeated. 
This procedure is continued until the temperature into the 
node is found to equal the saturation temperature or be 
above it. If the temperature is above the saturation temp- 
erature, then the m, and L parameters are changed to a value 
halfway between the current and previous values, and the 
process is continued until the temperature into the partial 
node exactly equals the saturation temperature. The values 


meet), T Sete and a are stored to be used later to 


cout fin 
develop weighted averages with the boiling portion parameters. 
The remainder of the node is solved in the boiling 


portion of the model. 


E. BOILING NODAL ANALYSIS 
Nodes in the boiling section are solved by the same 


iterative scheme developed in the superheating section. The 
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correlations utilized to determine the solution are differ- 
ent since flow in the boiling regime will involve two phases 
on the fluid-side. In two phase flow the inside heat trans- 
fer coefficient will change with quality, x, and the heat 
flux, q". The correlation selected for this model is the 
same as that used by Combs [Ref. 3] which was recommended 
by Tong [{Ref. 7] for both nucleate boiling and forced convec- 
tion. The analysis presented by Combs is presented here. 

The equation for the two-phase flow heat transfer coeffi- 


cient was given by Schrock and Grossman in Tong [Ref. 7] as 


TP af 
we = Bit + al) 
Me 15S tt 
The constants are given by Wright [Ref. 7]: B= BO 


A = Pace 107", ance = 0566.5 The heat transfer coefficient, 


assuming a totally liquid flow is 


k. &, CEG) Sa om 
anno ee ieee 
e oo Mo Q 


and the Martinelli parameter is 


The limits of the data for this correlation are: 
quality: ROIs On. 5 7 


meat flux: 6.0 i eee 10" BTU/hr-ft“, 
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In this model, heat flux in the boiling section will be at 
levels considerably below 6 «107 BTU/hr-£t* and the full 
range of quality must be considered. Therefore, in order to 
predict the performance of the Schrock and Grossman correla- 
tion at lower heat flux and quality above 0.57, a plot was 


made (Fig. 7), of the ratio /h, vs. quality where Ay is 


Dope 
the heat transfer coefficient from the Dittus-Boelter corre- 
lation for water at saturated liquid conditions. The corre- 
lation of Dengler and Addoms, as given in Tong [Ref. 7], 
indicates that low heat flux and high water-vapor mixture 
velocity favor the forced convection mechanism. After a 
review of the information available in Tong, it is not 
entirely clear whether the forced convection or the nucleate 
boiling mechanism is dominant for the relatively low flow 
velocities and low heat flux of this model. In any case, 
the Schrock and Grossman correlation applied at low heat 
flux represents a type of compromise between the entirely 
forced convection assumptions of Dengler and Addoms and the 
"mixed" assumption of Schrock and Grossman. 

A simplification assumed in the model is that the 
quality in the node is the average quality for the entire 
length of the node. Also, to expedite the numerical proce- 
dure, it is assumed that the node effectiveness is 0.25, 
initially. The effectiveness rarely attains a value higher 
than this for this type heat exchanger, so the assumption is 
valid. The value of the effectiveness is adjusted to the 


"true" value in the iterative process. Ina process involving 


Be) 





Meaiange Gt phase Ehe fluid is isothermal. In actuality, 
there is a slight temperature change through the boiling 
regime but for purposes of analysis this change is neglected. 
For this reason the fluid bulk temperature in any node in 

the boiling portion is equal to the fluid saturation tempera- 
ture. Therefore, to determine the amount of energy absorbed 
by the fluid in a node, the change in fluid enthalpy between 
the inlet and outlet is calculated. 

In the first node in the boiler, which is the node imme- 
diately before the superheater or a partial node if the node 
was split, the gas bulk temperature is again set equal to 
the gas inlet temperature. The gas-side heat transfer coeffi- 
cient is calculated using the same procedure as delineated 
in the superheater nodal analysis. Next, a Q for the node 


is calculated from 


rr 


Ceca © Mg Cae Ma Ag eelen 


where the effectiveness, ¢€, 1S initially assumed to be 0.25 


and Te sn ieewecual towtne flurd saturation temperature. A 


fluid inlet enthalpy is derived from 


Rey. = Ah 


= ag ie SSE (Q/m¢) 


where he saty = enthalpy of the fluid at saturated vapor. 


The fluid bulk enthalpy is 
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The average quality for the node is calculated from 


D eB 7 Be cael 


ay fsatv Ne satl 


where h 1s the enthalpy of the fluid at a saturated 


f satl 


liguid condition. The Martinelli parameter is determined 


from 


1 _; Xavg ‘a i ee Ob 
et 19% sug Ry nae 


where P,/0., are the ligquid/vapor densities at the saturation 
temperature and HL /U 9 are the liquid/vapor viscosities at 

the saturation temperature. Now the heat transfer coefficient 
Peeecocally liquid flow is calculated 


Os 


K d.G(l-x yOu Gale 
= } 


u 
Bi, =e! x 25 SE (2 
i ee Q 


From this the heat transfer coefficient is calculated for 


two-phase flow 





Ie 
h = hh, 6a ee 
Tet Me aay 3 are 
where h =h - h and the constants A, B andn 


1s | f satv f satl 


have the aforementioned values and G is the mass flow rate 
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per unit area. Finally, the node overall heat transfer 


Boerficient is calculated from 


U - ee pe 

OL le A. nid da.) A. i 
, _1P oo 1 6, _ip 1 

Dope 2m K Bh as L oy de Ny 








Whenever two phase flow is involved, the heat capacity of 
the fluid is assumed infinite. For this reason, within 
the boiling regime, Cin will always equal the heat capacity 


of the hot side, so, 


c 
min Gama 


and the Number of Transfer Units, NTU, is 


Un 4 Re 
NUM Tce 
min 
The node effectiveness is 
ce = 1 - exp(-NTU) 


Now, a revised Q for the node is calculated from 
Om = peeGc...._(<L oe ee | 


gluing | oi aaliee | Ele 


From this Q a revised gas outlet temperature is calculated 
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Ty out Ty in  '8/Cnin! 
With this gas temperature out, a new film temperature is 
calculated in the same manner as for the superheating section. 

The gas outlet temperature is stored and the film temp- 
erature is used to repeat the process for the node. At the 
end of the second time through, the gas outlet temperature 
derived is compared with the result from the first. If 
equality exists, a solution for the node has been found and 
the model proceeds to the next node. If inequality exists, 
the result of the current iteration is stored for comparison 
with the results of the following iteration. 

When a final solution is attained, the enthalpy of the 
fluid at the inlet to the node is checked against the 
enthalpy of saturated water at the heat exchanger pressure. 
If the inlet enthalpy is greater than or equal to the enthalpy 
of saturated water, the model proceeds to the next node and 
begins the solution process again. If the inlet enthalpy 
is less than the enthalpy of saturated water then the model 
proceeds to a node splitting routine similar to that used 
when transitting from the superheating regime to the boiling 
regime. 

The node splitting is again accomplished by halving the 
mass flow rate on the gas-side and halving the physical 
length, L, of the node. The remaining geometric parameters 


are recomputed and the solution process is reinitiated. The 
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mass and length adjustment procedure is continued until the 
inlet enthalpy is equal to the enthalpy of saturated water. 


When this is achieved, the values of T U_. and Q are 


Grout aa C1 
stored to be recombined as weighted averages with the heating 
portion of the node to give the overall nodal values. 


The remainder of the node is solved in the heating portion 


of the model. 


FF. HEATING NODAL ANALYSIS 

The procedure for analyzing a node in the heating portion 
of the model is identical to that used in the superheating 
portion. Each node is analyzed in an iterative process, and 
it makes no difference if the first node was a split node. 
The first node in the heating section will have a gas inlet 
temperature equal to the gas outlet temperature of the node 
directly below. The fluid outlet temperature will equal the 
fluid saturation temperature for the specified pressure. 

An outline of the calculations is as follows: 

1. Calculate the gas-side heat transfer coefficient 
using the node gas inlet temperature as the gas bulk 
mempenature:. 

2. Calculate the fluid-side heat transfer coefficient 
initially using the fluid saturation temperature as 
mie fluid bulk temperature. 

3. Calculate the overall heat transfer coefficient using 


the effectiveness-NTU method. 
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4. Determine the node fluid inlet temperature from 


- Zz = a Se a in | Me Cpt De Out 
ie ise! enc -m-C 


min i emong 


Seeeecalculate a first “guess" for Q for the node from 


Q = mM ¢ fogs out £ in) 


6. Calculate the node gas outlet temperature from 


Tgout ~ Tin ~ 7M, Sp) 


7. Calculate ue eee = ba + Oe ae) 2 


8. Calculate a node average film temperature from 


T = (Top tl le 
9. store the result of the first pass through the model 
and iterate until procedure converges to solution, 
utilizing gas outlet temperature for comparison at 
each pass. 
All properties of the gas and fluid are evaluated at the 
respective bulk temperatures for the node. The process 1s 
repeated node by node until the node at the fluid inlet to 


the heat exchanger is completed. The final node is either 
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node (1,1) for an odd number of passes, or node (1,10) for 


an even number of passes in the heat exchanger. 


G. HEAT EXCHANGER SOLUTION PROCEDURE 

When the procedure of the heating section has arrived at 
a solution for the final node of the heat exchanger, the 
fluid inlet temperature for the node is compared to the user 
specified feedwater inlet temperature. If equality exists, 
the model proceeds to an output routine where the matrices 
of water/steam temperatures, exhaust gas temperatures, 
overall heat transfer coefficients and individual heat 
thanster per node are printed. Additionally, the fluid-side 
mass flow rate, minimum pinch point temperature difference, 
total 9 for the heat exchanger and isentropic steam turbine 
power are calculated and printed. 

The pinch point temperature difference is defined in 
this model as the minimum difference in the heat exchanger 
between the gas inlet temperature and the fluid inlet temp- 
erature to any node. Consideration of the effect of a pinch 
point below approximately 25 deg. F is cause for concern 
only if the low pinch point is coupled with a low gas temp- 
erature. This combination is non-prohibitive to hot side 
corrosion. If the final node fluid inlet temperature is not 
equal to the user specified feedwater inlet temperature, then 
the fluid mass flow rate that was "guessed" at the beginning 


is adjusted and the entire model is run again. If the solution 
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temperature is too high the "guessed" fluid mass flow rate 

is lowered. With less mass going through the heat exchanger 
the flow velocity is decreased. A lower flow velocity through 
the heat exchanger enables more energy from the gas-side to 

be transferred to the fluid-side. A larger Q will increase 
the temperature differential between the fluid inlet and 
outlet which will yield a lower solution fluid inlet temp- 
erature the second time through the model. The process is 
Opposite, i.e., mass is added, if the fluid inlet temperature 
to the first node is found to be too low. 

The model will continue to adjust the mass after each 
iteration in smaller amounts until an exact match is made to 
the user specified feedwater inlet temperature. When equality 
is achieved, the model proceeds to the output routine where 
the above mentioned parameters are printed. 

Included, as Appendix A, are sample outputs of the model 
Meme ne uniform distribution, distribution 1 and distribu- 
tion 5 for all four of the selected gas mass flow rates. 

The run numbers at the top of the first page of each output 
were utilized for accounting purposes. The uniform distri- 
bution was used for runs 1, 10, 19 and 28; distribution 1 was 
used for runs 2, 1l, 20 and 29; and distribution 5 was used 
Gemerins 6, 15, 24 and 33. The first page of the output shows 
the inter-active questions asked by the model and the user 
responses. The following pages show the matrices developed 


by the model in the solution process. Comparison of the 
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matrices for the uniform and non-uniform distributions 
demonstrates the serious impact of the non-uniformity on the 
performance. The non-uniform distribution profiles can be 

seen as projected profiles on the results of each of the 
matrices. The final page of the output is a listing of 
additional parameters which are calculated by the model includ- 
ing the fluid mass flow rate, Q based on the individual 
Matrices, location of the beginning of boiling and super- 


heating within the WHRU and the isentropic horsepower. 
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feo CODSLON OF RESULTS 


ie Mier tAL TESTING 

The foregoing model was coded in FORTRAN IV and applied 
in a computer simulation on an IBM-3033 computer. 

To properly test the model and develop an understanding 
of the operation of the WHRU several runs were made utilizing 
a uniform gas distribution with a total mass flow rate and 
temperature corresponding to a gas turbine horsepower rating 
for approximately a 20 kt. speed. The conditions presupposed 
usage of a General Electric LM2500 Gas Turbine Engine, and 
gas mass and temperature data was taken from [Ref. 8]. Initially, 
the model was run with attention to pinch point temperatures. 
The number of passes was varied to find the largest heat ex- 
Changer which could be used and still maintain a pinch point 
in excess of 25 degrees fahrenheit. Under conditions of non- 
MirGerm f£low distributions the model definition of pinch point 
1s misleading since small gas mass flow rates through a 
particular gas path may yield a nearly zero pinch point when 
the overall performance of the heat exchanger is not affected. 
Inspection of the model output for the non-uniform distribu- 
tions will demonstrate this result. The preliminary runs 
allowed selection of a 23 pass heat exchanger with a 29.4 deg. 
F pinch point. Therefore, the 23 pass heat exchanger became 
the "design." The performance data tabulated for the LM2500 


engine [Ref. 8] was utilized to select four gas mass flow 
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rates and inlet temperatures which correspond closely to 
data points on a cubic power curve. The four combinations 


of mass and temperature chosen were: 


m. = 121.8 lbm/s, To in 7 844° F, SHP = 17500, NPT = 3000 rpm 
mM. = 99.1 lbm/s, To in = 755° F, SHP = 10000, NPT = 2400 rpm 
Se lbm/s, T, 5, = 702° F, SHP = 5000, NPT = 1800 rpm 
f= 00.7 lbm/s, T_.. = 694° F, SHP = 3000, NPT = 1200 rpm 

g G in 


Each of these was run in the model as uniform and various 
non-uniform distributions. 

To develop the non-uniform distributions, a normalized 
flow distribution table was constructed and is presented as 
Table I. This table shows the fraction of the total mass 
flow rate which is passing through each individual gas path 
for each of nine different flow distributions. The flow 
distributions were arbitrarily derived in an attempt to 
develop a spectrum of distributions that may enter the heat 
exchanger. A graphical display of each of the flow distribu- 
tions is shown in Figures 8A-8I. The total mass flow rate 
used was multiplied by each entry in the normalized flow 
distribution table to arrive at the mass flow rate through 
each of the ten gas paths for any of the nine distributions. 
This procedure allowed a comparison of the effects of non- 
uniform distributions on the performance of the heat ex- 


changer. All nine distributions were run for each of the 
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four chosen mass and temperature combinations. The actual 
mass flow distributions used are shown in Tables II-A--II-D. 
The combination of nine distributions and four different 
total gas flow rates yielded 36 runs. For the runs, the 


remaining parameters were: 


a) feedwater inlet temperature = 200 deg. F 
b) desired steam temperature = 650 deg. F 
c) steam operating pressure = 400 psia 
d) minimum allowed gas outlet temperature = 300 deg. F 
e) node length = d0 ft. 

f) number of tubes per row = 32 
fmsgeometric scaling factor = 1.0 


The model calculates the isentropic horsepower that would 
be developed by an ideal steam turbine. This power can also 
be viewed as the maximum power available from the developed 
Superheated steam exhausting to a fixed condenser pressure. 
Beeemetnod of calculation is to determine the enthalpy of the 
steam at the outlet pressure and temperature and then sub- 
tract the enthalpy of the steam after isentropic expansion 
to a condenser pressure of 2 psia. The difference in entnalpy 
was multiplied by the fluid mass flow rate, with the appro- 
priate conversions, to arrive at a horsepower value. The 
horsepower calculated for each run is tabulated in Table III. 
The table lists the results of each run of the model for all 
36 runs for each of 5 different steam outlet temperatures. 


Table IV is a similar table showing the solution values of 
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the fluid mass flow rates derived in each run. The two 

tables demonstrate the variation of fluid mass flow rate and 
power which is attained by varying the desired steam outlet 
temperature; and the effect of the non-uniform distributions 


on performance. 


Beer reCr OF NON-UNIFORM DISTRIBUTIONS 

The data from column 1 of Table III was plotted and is 
shown in Figures 9A-9E as power vs. gas flow rate. Figure 
9A shows the comparison of the uniform distribution with 
distributions 1 and 2. Figure 9B shows the comparison of 
Baemunttorm distribution with distributions 3 and 4, etc. 
From the graphs it can be seen that the distribution inflicting 
the most serious impact on the heat exchanger performance is 
Seeertoution 1 (see Figure 8-B). Distribution 1 is highly 
center concentrated and utilizes the area available for heat 
transfer the least efficiently. For comparison, distribution 
3 is a nearly uniform distribution and itinflicts a relatively 
Minor impact on the performance. Distributions 2 and 4 are 
ranked next in order of impact. These are also distributions 
where the mass is centrally concentrated, but to a lesser 
@eqgree than distribution 1. Distributions 5 and 6 are mirror 
images of one another and also have a zero mass flow rate 
through two of the gas paths. A zero mass flow rate indicates 
that the flow through that particular pass has no velocity. 


In a real Situation there will be a transient heat transfer 
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from the still mass of gas surrounding the tubes by natural 
convection until the gas temperature reaches the tube wall 

or fluid temperature. In any case, the amount of heat 
transferred in such a node is negligible compared to the 

nodes experiencing a gas flow rate. For the model, a zero 

gas flow rate in a node will be treated as having zero heat 
transfer, 1.e., a "dead node," therefore, the Q for the node 
was set to zero, the fluid inlet temperature and fluid outlet 
Gemperature were set equal and the gas outlet temperature 

was set equal to the gas inlet temperature. In reality, the 
gas temperature ina "dead node" will attain the temperature 
Of the fluid in the node but for pinch point calculation pur- 
poses the original gas inlet temperature was maintained. 

Also, since the model assumes no cross-mixing of the gas 
between gas paths as the distribution passes through the 

heat exchanger, this assumption ignores the potential tendency 
of the gas mass distribution to become uniform during passage. 
Distributions 5 and 6 are skewed to the left and right, 
respectively. The combination of the null gas flow paths and 
the skewedness of these distributions caused the second most 
detrimental effect on the heat exchanger. It should be noted, 
however, that whether the distribution was skewed to the right 
or left side had relatively no effect; that is, the results 
PerpeGrstribution 5 are nearly identical to those for distri- 
bution 6. The same result was true for distributions 7 and 


8, which are also mirror images of each other with left and 
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right skewedness. The effect of distributions 7 and 8 were 
very nearly the same as those for distribution 4. Figure 9E 
is a composite which demonstrates the effect of all the dis- 
tributions. The ranking as to effect yields: 

1. the uniform distribution gives the best performance 

Peeeaistr. 3, nearly uniform 

3. distr. 4, least center concentrated, nearly same as 

distr.'s 7 and 8 

4. distr. 2, intermediately center concentrated 

5. distr.'s 5 and 6, skewed with null flow paths 

6. distr. 1, highly center concentrated 

The ranking can be viewed as a measure of the efficiency 
of the distribution in utilizing the area available in the 
heat exchanger for heat transfer. As the amount of distor- 
tion of the flow increases, the quantity of energy trans- 
ferred decreases. 

Figures LOA-1lOE are graphical representations of the 
solution fluid mass flow rates vs. the gas mass flow rates 
Bemeeeae. Of the distributions; exactly as was done to develop 
Figures 9A-9E. As was to be expected, the results for Figures 
1OA-10E are identical to those for power vs. gas flow rate, 
Since power is directly proportional to the fluid mass flow 
rate. 

After these preliminary results, it was decided to con- 
duct the remaining runs utilizing only the uniform distribu- 


tions and the worst of the non-uniform distributions; l.e., 
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Gestrivoution 1. It was assumed that results for the remain- 


ing distributions would fall between the values of these two. 


Sore ocr OF HEAT EXCHANGER SIZE 

Next, an analysis was done to determine the effect on 
output power and fluid mass flow rate by varying the number 
of heat exchanger passes. Using the four chosen gas mass 
flow rates, runs were conducted for heat exchangers of 20, 
21, 23 and 25 passes. All runs were for 650 deg. F outlet 
steam at 400 pSia operating pressure. Therefore, the runs 
for the 23 pass heat exchanger were the "design" runs and 
the runs for 20, 21 and 25 passes were the "off-design" runs 
provided for comparison. The results of these runs are 
tabulated in Table V and graphically presented in Figures 
1l1A-11H as power and fluid mass flow rate, respectively, vs. 
the number of heat exchanger passes. Again, because of the 
relationship between power and fluid mass flow rate, the 
graphical results are similar in contour. The results show 
that as the size of the heat exchanger is increased, so will 
the power out. To provide the increased power the mass flow 
rate of the fluid must increase since for the different 
sized heat exchangers the inlet and outlet enthalpies of the 
fluid is not changing. The limit to the increase in size is 
the available energy in the gas. Increasing the heat ex- 
changer size to the point where sufficient energy is trans- 


ferred from the gas to the fluid to attain the minimum 
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acceptable exhaust gas temperature is the maximum size 
possible. 

The curves display an interesting result concerning the 
differences between the uniform distribution and the non- 
uniform distribution. As the size of the heat exchanger 
increases, the two curves are converging with the difference 
in power (or fluid mass flow rate) being less for the larger 
heat exchanger than for the smaller. From the data of Table 
V, it can be seen that, for the 121.8 lbm/s gas flow rate, 
the power from distribution 1 is 89% of the power attained 
from the uniform distribution for a 20 pass heat exchanger. 
This increases to approximately 91% for the 25 pass heat 
exchanger. The more rapid increase in power for the non- 
Uniform distribution is caused by the mass concentration in 
the center gas paths. For a particular size heat exchanger, 
the paths with a small gas mass flow will attain the temper- 
eeere or the fluid, usually in the boiling region. From that 
meemenon, these paths contribute little, or not at all, to 
the total heat transfer. The high concentration of mass 
in the center paths, however, contribute greatly to the total 
heat transfer because of the high velocity of the flow and 
the high heat transfer coefficients. These paths of high 
Goncentration will contribute until the gas exits the heat 
exchanger. By adding more passes to the heat exchanger 
it 1s as if the "dead areas" in the low mass flow paths were 


placed in the path of the high mass flows. The effect 1s a 
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significant increase in the heat transfer. The high heat 
transfer coefficients, mass velocity and temperature differ- 
entials all contribute to a greater encore nercase per 
pass than is achieved by the uniform distribution. As addi- 
tional passes are added, the performance of the non-uniform 
distribution improves, but it will never attain the level of 
the uniform distribution. For two heat exchangers of equal 
size, the uniform distribution gives the most efficient use 


of the available area for heat transfer. 


D. EFFECT OF STEAM OUTLET TEMPERATURE 

meen distribution, uniform and distributions 1-8, was 
Pum at each of the four gas flow rates for five specified 
steam outlet temperatures. The power and fluid mass flow 
Calculated by the model are presented in Tables III and IV. 
The results are also presented graphically in Figures 12-A 
eoeez—-6B for the uniform distribution and distribution l 
at each of the four gas flow rates. 

Both the tabular and graphical data shows a nearly uni- 
form increase in fluid mass flow rate with the linear de- 
crease in steam outlet temperature. The tabular data for 
power demonstrates a non-linear result, showing initially a 
decrease in power with a decrease in steam outlet temperature, 
and then increasing to attain a power higher than that for 
emeeu ceg. FF steam outlet temperature. 

The non-linearity of the power is caused by the non- 


linearity of the steam enthalpy as the steam outlet temperature 
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decreases. The combination of the change in enthalpy and 
the change in the fluid mass flow rate results in a mislead- 
ing impression of attaining a higher power at a lower steam 
temperature. The total power does increase but the specific 
power decreases. The increase in fluid mass flow between a 
650 deg. F steam outlet temperature and a 550 deg. F steam 
outlet temperature 1S approximately 7%. The increase in 
power at the lower steam outlet temperature is less than 0.2%. 
The penalty to utilize a lower steam temperature and main- 
tain an equivalent power output is a water/steam system which 
must be able to circulate 7% more fluid in the same time 
period. This result may have serious impact in the design 
when condenser size, circulating pump size, valves, piping 
and the steam turbine are taken into consideration. The 
steam turbine must also be considered with respect to the 
lower steam temperature having a lower quality through the 
latter stages of the expansion. The effects of overly "wet" 
steam can be disastrous to many turbines. The end result 

is that the sacrifice in power at the higher steam tempera- 
ture is justified by the decrease in fluid flow required 


and higher end point quality. 


PeeeeerreCr OF HEAT EXCHANGER PRESSURE 

ioe Uniform distribution and distribution 1 were run for 
each of the four gas flow rates with three different speci- 
fied heat exchanger operating pressures. The pressures were 


the "design" pressure of 400 psia and then 500 psia and 600 
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psia. The tabulated results of these runs are presented in 
Table VI, and graphically presented in Figures 13-A and 13-B 
for power and fluid mass flow rate, respectively, vs. heat 
exchanger pressure. Again, because of the direct correla- 
tion between fluid mass flow rate and power, the contours of 
the graphical results are similar. 

The curves show a decrease in power and fluid mass flow 
rate as heat exchanger pressure is increased. As the pres- 
Sure is increased by 50% above the design level, there is a 
4% drop in the system fluid mass flow rate and a 7.5% less 
in power. Clearly, the benefits derived from the decrease 
in system fluid mass flow rate does not justify the sacri- 
fice in system power. With respect to the analysis of the 
previous section, a compromise must be reached to achieve 
the maximum power from the system within the confines of a 
maximum acceptable fluid mass flow rate. 

The graphical results also show the effect of the non- 
Morrorm distribution for each of the four gas flow rates. 

As the gas flow rate decreases, the penalty for non-uniformity 
has a lesser effect. The percentage loss for the four flow 
rates are approximately 9, 8, 7 and 6 percent, respectively, 
from highest flow rate to lowest. The loss is approaching 
zero as the flow rate approaches zero. The loss in power is 
directly related to the decrease in fluid mass flow rate 
determined in the solution; and the decrease in fluid flow 
rate is a result of the ineffectiveness of the non-uniformity 


On the operation of the heat exchanger. 
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F, EFFECT OF GAS FLOW RATE ON HEAT EXCHANGER PERFORMANCE 

The data from the matrix of water/steam temperatures for 
the four gas flow rates as uniform distributions was plotted 
@ema cunction of tube length. The curves are shown in 
Figures 14-A--14-D. Figure 14-E is a plot of the water/ 
steam temperatures as a function of tube length for distri- 
bution 1 with a gas mass flow rate of 121.8 lbm/s. Figure 
14-F is a superposition of Figure 14-E on Figure 14-A. 

Analysis of Figures 14-A through 14-D shows the effect 
of the decrease in gas flow rate on the rate of heat trans- 
fer as the fluid passes through the heat exchanger. The 
higher the gas flow, the steeper is the slope of the curve 
in the superheat region. A steeper slope is indicative of 
a shorter tube length required to attain the desired super- 
heat temperature. THe effect, as the gas flow rate is 
decreased, is that more of the heat exchanger is dedicated 
to superheating. The length of the tube required to raise 
the water from a saturated liquid to a saturated vapor 
condition is nearly the same for all gas flow rates with 
Only a “left shift” on the plots caused by the superheat 
conditions. The slope of the curves in the heating region 
is the opposite of that of the superheat. As the gas flow 
rate decreases the slope of the curve in the heating region 
increases. 

Since, for all conditions, the feedwater inlet tempera- 


ture and steam outlet temperature are the same, all tubes 
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are continuous and of the same length, the difference in the 
Curves is attributable to the rate of heat transfer and the 


total amount of energy transferred. The rate ‘equation is 


al caus oe ae Te in? 


and this must equal the energy balance equation 
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The change in gas mass flow rate will obviously decrease the 
total amount of energy transferred since less gas in means 
less energy into the heat exchanger. If the curves were 
plotted in a third dimension with Q, the smallest Q would 
occur with the smallest gas flow rate. The forced equality 
with the rate equation means that, as the gas flow rate 
decreases, the rate of heat transfer must decrease. The 
figures show the total effect on the profile operation of 

the heat exchanger is of little consequence. The similarity 
of the curves is an indication of the "physical operation" 

of the heat exchanger under varying conditions. For the test 
runs, a decrease in gas flow rate was coupled with a decrease 
in gas inlet temperature. The two parameters changing and 
the similarity of the curves emphasizes the physical constancy 


of the design. 
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Pegs L4=— as the plot for distribution 1. The "ripple" 
effect in the heating and superheating regions is a result | 
of the non-uniformity of the gas mass distribution, but the 
superposition shown in Figure 14-F demonstrates the oblivious- 
ness of the physical properties of the heat exchanger to a 
Mereration Of input. 

G. EFFECT OF NON-UNIFORM MASS DISTRIBUTION ON EXHAUST GAS 

TEMPERATURE 

For all runs of the model the temperature distribution 
of the exhaust gas entering the heat exchanger was assumed 
uniform. To determine the effects of a non-uniform gas mass 
distribution on the exhaust gas temperature profile, a plot 
was made of the gas temperatures across several of the heat 
exchanger passes. The matrix of exhaust gas temperatures 
was used and data was plotted for both a uniform distribu- 
tion and distribution l. The gas mass flow rate was 121.8 
lbm/s. Figures 15-A and 15-B show the results. 

Figure 15-A demonstrates that with a uniform mass and 
temperature distribution at the inlet to the heat exchanger, 
the distribution is virtually undistorted by passage through 
the WHRU. This result must be considered ideal with respect 
to the assumption of no cross-mixing of the gas during passage. 
There is a smooth transition of energy from inlet to outlet 
which accounts for the best performance of the heat exchanger 
under this condition. The effect of the non-uniformity 1s 


clearly seen in Figure 15-B. The temperature distribution is 
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rapidly distorted to become almost a replica of the gas flow 
distribution. This is a result of the disparity of heat 
transferred and heat transfer rates from node to node in 

the WHRU. A slight smoothing effect is seen in the boiling 
portion, but distortion is again accentuated in the heating 
section. In areal situation there will be some non- 
uniformity of temperature as it enters the heat exchanger. 
The analysis shows that the condition will be worsened to 
the degree of the mass non-uniformity with the subsequent 
loss of performance in the WHRU. 

The distortion of the exhaust gas temperature distribu- 
tion forces consideration of pinch point. When the tempera- 
ture differential between the exhaust gas and the tube 
becomes less than approximately 25 deg. F coupled with a gas 
temperature less than approximately 300 deg. F, the WHRU may 
experience localized hot-side corrosion due to sulphur products 
in the exhaust gas. It should be noted that this condition 
Will exist only near the top of the heat exchanger near the 
exhaust gas exit. If the degree of flow distortion is 
severe and maintained for a long period, the effects can be 


disastrous. 
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IV. CONCLUSIONS 


The waste heat recovery unit iS a practical method of 
recouping thermal energy from gas turbine engine exhaust 
gases for fuel economy improvement. In the design and 
installation of a RACER system several factors must be taken 
into consideration. 

1. Mass Flow Distribution. The performance of a waste 
heat recovery unit is adversely affected by a non-uniform gas 
mass distribution; and the degradation is directly propor- 
tional to the degree of non-uniformity. To ensure acceptable 
results, the system should have flow distribution control 
devices installed between the gas turbine and the WHRU. 

2. Operating Pressure. The operating pressure of the 
system affects the power output inversely. A higher pressure 
System requires less fluid to deliver the same power but the 
percent reduction in fluid required does not justify the 
percent loss in power. Since the system is not a main pro- 
pulsion system requiring a large volume of fluid, the model 
Suggests using a low pressure system. 

3. Steam Outlet Temperature. Reducing the outlet temp- 
erature required forces an increase in fluid mass flow rate 
to maintain the same power out. High temperature steam, how- 
ever, has a greater wear effect on system components such as 
valves, seals and piping. The model suggests a steam tempera- 
Mane de High as practical should be used, within material 
constraints. 
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4, Heat Exchanger Size. Power and fluid mass flow rate 
are directly related to heat exchanger size. The size of 
the WHRU should be as large as the available ‘space and 
weight constraints allow, or as large as necessary to pro- 
duce the power required. Consideration must be given to the 
system as a whole since, as the heat exchanger becomes larger, 
so does the piping, valves, pumps and condenser necessary to 
Operate. 

5. Speed. Power developed is a function of the energy 
available in the exhaust gas flow. At higher gas turbine 
speeds more horsepower is available from the WHRU to assist 
the main propulsion system. Utilizing the latest NAVSEA 
speed profile, the heat exchanger should be sized in accordance 
with the ship speed most commonly used. 

All data suggests a WHRU as large as physically possi- 
ble, operated at a low steam pressure, high steam temperature, 
and with flow distribution control devices installed provides 


the best performance. 
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V. RECOMMENDATIONS FOR FURTHER RESEARCH 


A. WHRU MODEL IMPROVEMENT AND EXPANSION 

The control imposed on the WHRU model in this study was 
adjustment of the steam flow rate to maintain a user speci- 
fied steam outlet temperature. This forces the steam system 
to respond to a variation of gas inlet conditions by "throt- 
tling" the feedwater system to maintain the required steam 
eemperature. 

This is only one of several types of control systems 
which may be used in a RACER system installation. Two 
possible alternative systems are: one which maintains a 
constant steam flow rate while allowing the steam tempera- 
ture to fluctuate, and one in which pressure is allowed to 
"float" with steam demand. In this second type system, the 
flow resistance of the system is used, and as the fluid mass 
flow rate decreases (Steam demand decreases), the pressure 
required decreases. Knowing the inter-relationship will 
yield a "floating” out as well. Investigation into the 
performance of a WHRU with these types of control could lead 
to an improved output. The investigation of all types of 
control should consider the effects on the entire system 
because, the type of control may dictate the type steam 
turbine which must be used. 

To expand the understanding of the effects of non-uniform 


gas caistributions on the WHRU performance, the current 
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model could be enlarged to allow a variation of the input in 
an additional dimension. A two dimensional non-uniform input 
will yield a three dimensional picture of the temperature 
distributions as the gas passes through the heat exchanger. 
An analysis of this type would help design the flow distri- 
bution control devices to eliminate severe trouble spots 
contributing to hot side corrosion and to maximize the WHRU 
Pemrormance. 

In addition to the expansion to a three dimensional 
analysis, the current model could be improved by including 
the following items. 

1. Fluid-Side Pressure Drop 

The water/steam-side pressure drop may be a signifi- 
cant percentage of the system operating pressure. Calculation 
of the pressure drop and its effect on system performance 
Should be included. 

2. Additional Geometries 

Additional fin-tube configurations should be con- 
Sidered for possible improvement of the heat transfer charac- 
teristics and the weight and space requirements of the WHRU. 

eee rOulLiIng 

Consideration of inside and outside heat transfer 
surface fouling should be given when investigating possible 
fin-tube configurations. 

4. Alternate Working Fluids 

Working fluids, other than water, should be investi- 


gated for possible enhanced thermodynamic characteristics at 
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a minimum cost in system maintainability, while maintaining 
safety to operating personnel. 
a OlmumErormememperature Distributions 

The current model could be utilized to determine the 
effects of non-uniform temperature distributions. The gas 
is at a uniform temperature as it exits the gas turbine, but 
it probably is not uniform at the inlet to the WHRU. Inves- 
tigation into the effects of non-uniformity of both gas mass 
and temperature should be accomplished. 

The model should be incorporated into a complete 
RACER system analysis, including models of the steam turbine, 
condenser and circulating pump. Variations of gas distribu- 
tion input can then be analyzed to determine total system 
effect. Understanding the total system effect and the inter- 
relationship of the variables could then be used to optimize 


Bae System. 
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SFaleLi Sy 


GAS FLOW DISTRIBUTIONS (m, = 121.8 LBM/S) 


GAS MASS FLOW RATE = 121-8 LEMSS 
CaS TEMPERATURE AT HEAT‘EXCHANGER INLET = 88% DEG. F 
‘case 575 {REFS ) SHP 2= 175590 xNPT = 3000 
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TABLE II-B 
GAS FLOW DISTRIBUTIONS (m,, = 99,1 LBM/S) 


GaS “aSS FLOW RATE = 9922 Lanes 
GAS TEMPERATURE AT HEAT EXCHANGE? INLET = 755 DEGe F 


CASE S6t (REF 8) S¥P = 10000 NPT = 2400 
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TAS Gee LC 


GAS FLOW DISTRIBUTIONS (m1, = 77.4 LBM/S) 


GaS MASS FLOW RATE = 772% LBMSS 
GAS TEMPERATURE AT HEAT EXCHANGER INLET = 702 DEG. F 
cCaSE 553 (REF 8) SHP £ $000 NPT = 1800 
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TABLE II-D 


GAS FLOW DISTRIBUTIONS (mm, = 66.7 LBM/S) 


GAS MASS FLOW RATE = 67.7 LBMZS 


GAS TEMPERATURE AT HEAT EXCHANGER INLET = 69% DEG. F 
CASE 505 (REFS) SHP = 3009 NPT = 1200 
DIST. GAS PATH 
NO. 
1.15 re cir. 20 | 17.20 ee 4.50 Bie 32 
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TABEE ILI I 


ISENTROPIC POWER FOR VARIOUS GAS FLOW DISTRIBUTIONS, 
GAS FLOW RATES AND STEAM OUTLET TEMPERATURES 


DISTR. STEAM TEMPERATURE 
NO. 
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TABLE IV 


FLUID MASS FLOW RATE FOR VARIOUS GAS FLOW RATES, 
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TABLE Vs POWER AND FLUID FLOW RATE 
FOR VARIOUS SIZE HEAT EXCHANGERS 
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TABLE V1. PQWER AND FLUID FLOW RATE FOR VARIOUS 
STEAM DPERATING PRESSURES 
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LAYOUT OF MARINE GAS TURBINE WASTE-HEAT STEAM GENERATOR: ONCE-THROUGH, 
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Figure 1. Waste Heat Recovery Unit Layout 


76 





= Number of Heat Exchanger Passes 


General Node (m,n) 


WS. 


2,2 1.2) Em 


NI. 





aa A fo 
HW NNN Me 


Gas Paths 


Figure 2, Heat Exchanger Nodal Arrangement 
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Figure 6. Waste Heat Recovery Unit Tube Layout 
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APPENDIX A 


Sample runs of the computer output are provided for a 
comparison of the effects of the non-uniform gas mass 
distributions on the temperatures within the heat exchanger 
and the effect on overall performance. 

Runs 1, 2 and 6 utilized a total gas flow rate of 
ies lom/s with the uniform distribution, distribution 1 and 
Gistribution 5, respectively. Runs 10, 11 and 15 utilized 
the same distributions but with a total gas flow rate of 
99.1 lbm/s. Runs 19, 20 and 24; and runs 28, 29 and 33 
also utilized the same distributions but with mass flow 


rates of 77.4 lbm/s and 66.7 lbm/s, respectively. 
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